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Abstract 



Navier equations are used to describe the deformation of a homogeneous, isotropic and 
Unear elastic medium in the absence of body forces. Mathematicahy, the system is a natural 
vector (field) 0(n, M)-invariant generalization of the classical Laplace equation, which physically 
describes the vibration of a string. In this paper, we decompose the space of polynomial solutions 
of Navier equations into a direct sum of irreducible 0(n, M)-submodules and construct an explicit 
basis for each irreducible summand. Moreover, we explicitly solve the initial value problems for 
Navier equations and their wave-type extension — Lame equations by Fourier expansion and Xu's 
method of solving flag partial diff'erential equations. 



1 Introduction 

Classical Laplace equation 

'^xixi ~l~ '^X2X2 ~l~ ■ ■ ■ ~l~ Ux„x„ (1-1) 

is one of most fundamental partial differential equations in mathematics and physics, 
whose solutions are called harmonic functions. Physically it describes the vibration of a 
string. Its more general form on Riemmanian manifolds is the main object in harmonic 
analysis. A fundamental algebraic characteristic of the above equation is its invariance 
under the action of the orthogonal group 0{n, M), that is, the space of harmonic functions 
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forms an 0(n, R)-module. Denote by Ak the space of polynomials of degree k and by Tik 
the space of harmonic polynomials of degree k. The classical harmonic analysis says that 

Ak^nke{xl + xl + --- + xl)Ak-2 (1.2) 

and Tik forms an irreducible 0(n, IR)-submodule. By induction, the above conclusion 
gives a decomposition of the polynomial algebra into a direct sum of irreducible 0{n, M)- 
submodules. Another closely related fundamental equation is the wave equation: 

The solutions of the initial value problems for the Laplace equation (1.1) and the wave 
equation (1.3) are elementary known facts appeared in many textbooks of partial differ- 
ential equations. 

A partial differential equation of fl,ag type is a linear differential equation of the form: 

{dl + fld2 + /24 + ■ • ■ + fn-ldn){u) = 0, (1.4) 

where di,d2,---,dn are certain commuting locally nilpotent differential operators on the 
polynomial algebra ]R[xi, X2, Xn] and /i, fn-i are polynomials satisfying 

drifs) = if r > s. (1.5) 

Flag partial differential equations naturally appear in the problem of decomposing the 
polynomial algebra (symmetric tensor) over an irreducible module of a Lie algebra into 
the direct sum of its irreducible submodules. Many important linear partial differential 
equations in physics and geometry are also of flag type. Xu [XI] used the grading technique 
in algebra to develop methods of solving such equations. In particular, he found new 
special functions by which we are able to explicitly give the solutions of the initial value 
problems of a large family of constant-coefficient linear partial differential equations in 
terms of their coefficients. 

To avoid the confusion with time variable t, we use up-index "T"' to denote the trans- 
pose of a matrix (vector) throughout this paper. Navier equations 

iiA{u) + (ii + t2) ( • V) (u) = (1.6) 

are used to describe the deformation of a homogeneous, isotropic and linear elastic medium 
in the absence of body forces (e.g., cf. [G], [M], [T]), where u is an n-dimensional vector- 
valued function, A = dl_^ + d^^ -\ h d^^ is the Laplace operator, V = {8x^,8x2: dx„) 

is the gradient operator, ti and L2 are Lame constants with li > 0, 2ti + ^2 > and 
I'l + i'2 7^ 0. In fact, • V is the well-known Hessian operator. Mathematically, the 
above system is a natural vector (field) 0(n, ]R)-invariant generalization of the classical 
Laplace equation. Our first objective in this paper is to decompose the space of polynomial 
vectors into a direct sum of irreducible 0(n, M)-submodules in terms of the irreducible 



0(n, M) -sub modules included in the subspace of homogeneous polynomial solutions of 
Navier equations. Unlike the case of the classical Laplace equation, the subspaces of 
homogeneous polynomial solutions of Navier equations are not irreducible. 

The wave-type extensions of Navier equations are the well-known Lame equations 

Uu = ^^Mu) + {V^-V){u). (1.7) 

f'l + '-2 

Our second objective in this paper is to solve the initial value problems for Navier equa- 
tions and Lame equations by Xu's method of solving flag partial differential equations in 
[XI] and his matrix-differential-operator approach in [X2]. Below we give a more detailed 
technical introduction. 

Olver [O] (1984) found the symmetry group and conservation laws for Navier equations. 
Moreover, Ozer [O] (2003) obtained some exact solutions of Navier equations with n = 3 
by means of Lie point transformations. Rodionov [R] (2006) studied finite solutions of 
Lame equations by holomorphic expansions. From an algebraic point of view, it is more 
convenient to deal with the orthogonal Lie algebra o{n, M) than the orthogonal Lie group 
0{n, M). They are equivalent under the exponential map from the Lie algebra to the Lie 
group in our case. 

Similar results as (1.2) hold for some other Lie algebras, too. Denote by Er^s the square 
matrix with 1 as its (r, s)-entry and as the others. Recall the special linear Lie algebra 

n-l 

Sl{n, M) = ^ REr,s + J2 ^i^r,r - ^^r+l.r+l)- (1-8) 
r^s r=l 

Let A be the polynomial algebra in xi, yi, y„. Define a representation of sl{n, R) 
on A via 

Er,s\A = Xrd^^ - Vsdy^. (1.9) 

Denote by N the additive semigroup of nonnegative integers. Define 



Set 



and define 



= for a = («!,..., «„) e N". (1.10) 



= for£i,£2eN (1.11) 

a,/3eN"; \a\=h, \l3\=i2 



'^iiA = {ue Ai,,e2 I U^iyi + Ux2y2 + " " " + U^nVn = 0}- (1-12) 

Xu [XI] proved that Ti-i^/^ forms an irreducible sZ(n, ]R)-submodule and 

-^tiA = T^liA ® (^ll/l + ^2|/2 H h Xnyn)Al^-l,l^-l, (1.13) 



where Xiyi+X2y2 + ■ ■ • + Xnyn is an sl{n, R)-invariant. More importantly, an explicit basis 
for each Ti-i^^/^ constructed in [XI]. 



Now we assume that A is the polynomial algebra in Xi, ...,Xy. There exists an action 
of the simple Lie algebra Q'-''^ of type G2 on A, which keeps + 0:2X5 + x^xq + X4X7 
invariant (e.g., cf. [H] and [XI]). Again we denote by Ak the subspace of polynomials of 
degree k in A and define 



Xu [XI] proved that Hk forms an irreducible ^"^^-submodule and 

Ak = Hk® {xj + X2X5 + X3X6 + X4Xj)Ak-2- 



(1.14) 



(1.15) 



Moreover, an explicit basis for each Hk was constructed in [XI]. Furthermore, Luo [L] 
generalized the result (1.2) to certain noncanonical polynomial representations of o{n, R). 
To state our results in this paper, we denote 



X 



fix) 



( h{x) \ 

\ fn{x) j 



1.16) 



i = A with Ak = {/I e A}, (1.17) 

A;=0 

where Ak is the space of polynomials in xi, x„ of degree k. Moreover, we define 

Hk = {/e Ak I iiA(/) + (ti + i2)(V^ • V)(/) = 0}. (1.18) 



The action of the orthogonal group 0(n, M) on A is defined by 

T(/(f)) = {Tf){T-\x)). 



(1.19) 



Furthermore, we denote = (0, 0, 1, 0, 0)^ and h — {11+12)/ ii- Then / = Y^^=i fr^r 
We define linear maps ip,(pi,(f2 : M[xi, x„] — > .4 by 



;i.20) 



s=l 



(pi{xi,Xi^ ■ ■ -XiJ = ^Xii ■ ■ -Xi^.iXi^+i • ■ ■Xi^[r(^xl)<;i^-{2r + n-2)xi^^xi<;i] (1.21) 

s=l /=1 1=1 

and (/32 = (a^i H h xDijj. 

Again we let Hk be the space of harmonic polynomials of degree k. Set 



1=1 



Hk,! = V'(^ik+i), ^fe,2 = { I // e Hk, J^xifi = 0}, (1.22) 

(1.23) 



(2fe + n-2)(fc + n^3)(fc-l) 
^^'^ = I + 2(6-i(2)t + n-4) + )t-l) ^^j 



Denote dr,s = Xgdx^ — Xrdx^ and set 

/ ^3,4 ^4,2 C?2,3 \ 

C?2,4 C?4,l (ii,2 
V C?3,2 rfi,3 ^2,1 / 

With n = 4, we define 

nk,2±^{fenk,2 I p/ = ±(A; + i)/} 

Main Theorem 1. Le^ n > 3 be an integer. The subspace of polynomial solutions 
'Hk = Hk,! ® 'Hk,2 ® 'Hfe,3 anc? 

Ak = Hk ® {xj + ■ ■ ■ + xl)Ak-2- (1.24) 

Moreover, the subspaces ii.k,i, '^fe,2 (^t- 4) andTik^s are irreducible 0{n,M.)-submodules. 
When n — A, 7i.k,2 — 'Hk,2+ ® 'Hk,2- and 'Hk,2+, 'Hk,2- are irreducible 0{A, W)-submodules. 

Furthermore, an exphcit basis is constructed for each of the above irreducible sub- 
modules. Xu's Lemma in [XI] on the polynomial solutions of constant-coefficient partial 
differential equation of certain type is used to prove the equation Hk — 'Hk,i®'Hk,2®'Hk,3- 
Our above result is exactly a natural vector generalization of those scalar ones in (1.2), 
(1.13) and (1.15). Indeed, it does not only reveal the internal symmetry of Navier equa- 
tions but also gives the explicit realizations of certain irreducible 0{n, R)-modules. 

Our second main result is obtaining explicit exact solutions of Navier equations (1.6) 
subject to the initial conditions 

u{0,X2,...,Xn) ^ go{x2,...,Xn), Uxi{0, X2, Xn) ^ gi{x2, Xn) (1.25) 

for ], and explicit exact solutions of Lame equations (1.7) subject to the 

initial conditions 

{{(O, Xi, Xn) = hQ{xi, X„), ^^(0, Xi, Xn) = hi{xi, x„) (1.26) 

for Xr G [— 6r,&r], where gi, g2,hi,h2 are vector-valued continuous functions, and ar,br 
are positive real constants. Xu's method of solving the initial-value problems of flag 
partial differential equations in [XI] and his matrix-differential-operator approach [X2] 
play fundamental roles in obtaining our solutions. Moreover, Fourier expansions are used. 

Section 2 is devoted to the proof of Main Theorem 1. In Section 3, we construct 
an explicit basis for each irreducible submodule included in the solution space of Navier 
equations. Moreover, an uniform explicit basis of the homogeneous polynomial solutions 
of Navier equations is obtained by a method of Xu in [XI] , whose cardinality was pre-used 
in Section 2 to prove the completeness of polynomial solutions for a technical reason. In 
Section 4, we solve the above mentioned initial value problems. 



(1-24) 



(1-25) 



2 Polynomial Solutions and Representations 



In this section, we will study the homogeneous polynomial solutions of (1.6) in detail. 

As we mentioned in the introduction, studying 0{n, R)-representation structure of the 
polynomial solutions is equivalent to studying their o{n, ]R)-representation structure via 
the exponential map from the Lie algebra to the Lie group. Recall that Er^s is the square 
matrix with 1 as its (r, s)-entry and as the others. The orthogonal Lie algebra 

n 

o{n,R)^J2^(^r,s-Es,r). (2.1) 

r,s=l 

Its action on the space A of polynomial vectors (of. (1.17)) is given by 

{Er,s - Es,r){f) = {Xrd^^ - Xsd^^){f) + /^^^ - (2.2) 

where = (0, 0, 1, 0, 0)"^ as in the introduction. Moreover, the elements of o{n, M) 
map solutions of (1.6) into solutions. That is, 

^{Hk) C Hk for C G o(n, M), A; G N (2.3) 

(cf. (1.18)). Note that f+gi is a complex solution of Navier equations (1.6) if and only if / 
and g are real solutions. In order to use the representation theory of the Lie algebras over 
the complex field C, we need to complexify our vector spaces, indicated by the subindex 
of C. The complexification of o(n, R) is o(n, C). We extend the representation (2.2) of 
o(n, R) on A to that of o{n, C) on C-bihnearly. 

Denote by Mmxm(^) the algebra of m x m matrices with entries in the field F. Set 




o'(2m+ 1,C) = <( I a A Ai \ \ A e M^^^{Cy,Ai,A2 e o{m,Cy,a,b e 

(2.4) 



and 



o'(2m,C) = |(^ J^V ) l^eM^xm(C)Mi,A2eo(m,C)|. (2.5) 

Note that o'(3, C) is the standard form of complex simple Lie algebra of type Ai, o'(4, C) 
is the standard form of complex semisimple Lie algebra of type Ai (BAi and o'(6, C) is the 
standard form of complex simple Lie algebra of type A^ (or D^). Moreover, o'{2m + 1, C) 
with m > 2 is the standard form of the complex simple Lie algebra of type B.^ and 
o'(2m, C) with m > 4 is the standard form of the complex simple Lie algebra of type D^n. 
If n = 2m + 1 is odd, we take 

m 

= C(il^r+l,r+l ~ Em+r+l,m+r+l) (2.6) 



r=l 



as a Cartan subalgebra of o'{2m + 1, C), 



{Er+l,s+l — Efn^s+l,m+r+l, -E'r-+l,m+s+l ~ -£'s+l,m+r+l 

Er^^i - Ei^m+ri \ I < T < s < m; 1 < Ti < Ul} (2.7) 

as positive root vectors, and 

{Es+l,r+l ~ Em+r+l,rn+s+ly -E'm+s+l,r+l ~ -^m+r+l,s+l 

Em+ri,i - Ei^n \ 1 < r < s < m;l < ri < m] (2.8) 

as negative root vectors. 

For 1 < r < m, we define tlie linear function Sr on by 

^r{Es+l,s+l ~ Efn+s+l,m+s+l) — ^r,s- (2-9) 

Then H"^^^^ = span {ei, . . . , Sm} and $ = {±Sr, ±(£r ± £s) | 7^ s} forms tlie root system 
of o'{2m + 1, C). Denote by ar — Sr — £r+i for r = 1, . . . , m — 1 and am = £m- As we 
know, {«!, . . . , a„i} is the simple root system. Moreover, we define a symmetric bilinear 
form (•, •) on by 

{er,es) = Sr,s. (2.10) 
Furthermore, we denote by Ai, . . . , the fundamental dominant weights in H^^, i.e. 

[Xr, as) = r- = dr,s- 

In particular, 

£1 = Ai, £2 = — Ai + A2, . . . , £m-l = + Am-1, £m = —^m-l + 2Ato- (2.12) 

When n — 2m is even, we take 

m 

— C{Er,r — Em,+r,m+r) (2-13) 

as a Cartan subalgebra of o'(2m, C), 

{Er,s ~ Em+s,m+ri Er,m+s ~ Eg^m+r | 1 < T < S < m} (2-14) 

as positive root vectors, and 

{Es,r ~ Em+r,m+s, Em+s,r ~ Em+r,s | 1 < T < S < m} (2-15) 

as negative root vectors. Moreover, we define Sj. e H*j-,^ by 

£r{Eg g Em,-\-s,m+s} ~ ^r,s- (2-16) 



Similarly, H"^^ = spanjei, . . . ,em} and $ = {±(£^±£5) \ r ^ s} forms the root system of 
o'(2m, C). Moreover, we define a symmetric bilinear form (•, ■) on by (2.10). Denote 
by ar = Sr — Sr+i for r = 1, . . . , m — 1 and am = Sm-i + £m- Then {ai, . . . , is the 
simple root system. Moreover, Ai, . . . , A^^ are the fundamental dominant weights in H}j^ 
defined by (2.11). In this case, 



£1 — Ai, £2 — —^1 + ^2, ■ ■ ■ , £m-2 — + 

Denote by 1^ the m x m identity matrix. Set 

1 + i 



if n = 2m + 1, and 



K 



1-i 



^-^m -^m 



K 



\ — i / iT T 
I il 



when n — 2m. Then X is a symmetric matrix, 

/1-i 



1+i 





V 



-il. 



if n = 2m + 1, and 



K 



-1 



1 + i 



m 

Im Urn 



-il I 



when n — 2m. Furthermore, we have a Lie algebra isomorphism a : o'{n,C) 
given by 

a(X) = K-^XK for X e o'(n, C). 
We define a representation p of o'(n, C) on A by 

p(X)(/) = <j(X)(/) for X e o'(n,C). 

So A forms an o'{n, C)-module. 

In order to study the o'{n, C)-module structure of ^c, we let 



\yn J 



1+i 



(2.17) 
(2.18) 



(2.19) 



(2.20) 



(2.21) 



(2.22) 
o(n, C) 

(2.23) 
(2.24) 



Define two representations pi and p2 of g'/(n, C) on Ac = C[a;i, 



(2.25) 
by 



(2.26) 



Writing K = {br,s)nxn and K ^ = (c^,s)„xn, we have 



and 



Thus 



K Ep qK ^ ^ Cs^pbq^rEg^r 
r,s=l 



n n 



dxs dyg 



r,s=l 



p^{X) = p2{K-^XK) for X e gl{n, C). 
RecaU = (0, 0, 1, 0, 0)^. Furthermore, we set 



/ 








\ 










2 






\ 











Then 



C ^s- 



s=\ 



We define two representations p\ and p2 of C) on Ac by 



P2(X)(^gr<;r) ^^P2{X)(gr)<;r + (gi, ...,gn)X^ 



r=l 



r=l 



for fr,gr e and X e gl{n,C). By (2.27)-(2.31), we have 

= p2{K-^XK) for X e ^/(n, C). 
According to (2.2), (2.23), (2.24) and (2.34), we obtain: 
Lemma 2.1 For X e o'(n, C), p{X) = pi(X). 

Denote by Ac,k the subspace of the polynomials with degree k in ^c- Set 

n n n 



V 



r=l 



s=l 



= {/ e ^cfe I (a,\ + 2 ^ a,,^,a,^^,^j(/) = 0} 



if n = 2m + 1, and 



K = {/ e ^Cfc I (E = 0} 



if n = 2m. It can be verified that 

n 

m^C^^Hk^ife Ac,k I iY. = 0} (2-38) 

r=l 

(also see (2.68)). Now V forms an o(n, M)-submodule of A that gives the basic rep- 
resentation with the canonical basis {<?i, <^n} (cf. (2.2)). Moreover, Vc forms an 
o'{n, C)-submodule of Ac that gives the basic representation with the canonical basis 
{Ki,...,K,n} (cf. (2.30)). Furthermore, A forms an o(n, M)-module with the represen- 
tation p2\o(n,R), and forms an o'(n, C)-module with the representation pi|o'(n,c) (cf- 
(2.26)). As a real o(n, M)-module, Ak = AkV ^Ak^wV (cf. (2.2)). On the other hand, 
Ac,k — AckVc — Ac,k 'S>c Vc as a complex o'(n, C)-module by the above lemma. We want 
to decompose Ak as a direct sum of real irreducible o{n, M)-submodules via decomposing 
Ac,k as a direct sum of complex irreducible o'(n, C)-submodules. To study the complex 
o'(n, C)-module ^c,fc0c^c, we recall a fact about the tensor product of modules for a 
finite-dimensional complex semisimple Lie algebra Q. 

Denote by V{X) a finite-dimensional irreducible ^-module with highest weight A, by 
A the weight lattice of G and by A+ the set of dominant weights. Moreover, V^{X) stands 
for the weight space oiV{X) with the weight u. 

Lemma 2.2 (e.g., cf. [GE]) The irreducible representations occurring in V{X) ^V{ii) 
have highest weights of the form fi + u E A+, where u is a weight in V{X). Moreover, the 
multiplicity ofV{ii + u) in V{X) V{ii) is less than or equal to dim T4(A). 

The weights of o'(2m, C)-module Vc = ^(Ai) with m > 3 are 

ei^ e2> £m-2 y- £m-i ^ ^ >- -£m-i > ^ "^i- (2-39) 

Moreover, the weights of o'(2m -|- 1, C)-module Vc — ^('^i) with m > 2 are 

Siy €2)- ... em-2 >- £m-l £m > Sm > £m-l '>-■■■> £2 > £l- (2.40) 

The weights of o'(4, C)-module Vq = V{Xi + As) are 

A1 + A2 ^7^^ + ^^ y-X,-X2. (2.41) 

The weights of o'(3, C)-module Vq = 1^(2Ai) are 

2Ai ^ ^ -2Ai. (2.42) 

Each weight in Vc is with multiplicity 1 (cf. (2.9) and (2.16)). 
It is well known that 

2m 

Ac,k ^n'k®{yl + 2j2 yr+iym+r+i)Ac,k-2 (2.43) 

r=l 



if n = 2m + 1, and 



2m 



Ac,k = K ® (Yl yrym+r)Ac,k-2 (2.44) 
r-=l 

if n = 2m. Moreover, yi+2 X]r=i Z/r+iZ/m+r+i is an o'(2m+l, C)-invariant and J2l=i VrVm+r 
is an o'(2m, C)-invariant. Furthermore, H'^ = T^(A;Ai) if n > 4, H'^ ^ y(2A;Ai) if n = 3 
and n'j, = y (A;Ai + k\2) if n = 4. 

Lemma 2.3 As o'{n, C)-modules with n>7, 

n'kVc ^ H'k ®c Vc ^ V{{k + l)Ai) © V{{k - l)Ai + A2) © V{{k - l)Ai). (2.45) 

a) For o'(6,C), 

= K ®c = ^((fc + l)Ai) © V{{k - l)Ai + A2 + A3) © V{{k - l)Ai). (2.46) 
in) For o'(5, C), 

n'kVc ^ n'k ®c Vc ^ v{{k + i)Ai) © v{{k - i)m + 2A2) © v{{k - i)Ai). (2.47) 

iv) Foro'{4,C), 

KVc = K Vc 

^ V{{k + l)(Ai + A2)) © V{{k - l)Ai + {k + 1)A2) © V{{k + l)Ai + {k- 1)A2) 
©y((A;-l)(Ai + A2)). (2.48) 

v) For o'(3,C), 

n'kVc = Tifc (8)c V^c = V^(2(A; + l)Ai) © V^(2A;Ai) © V{2{k - l)Ai). (2.49) 

Proof, i) First we consider the case that n — 2m is even with m > 4. Lemma 2.2, (2.17), 
(2.18) and (2.39) tell us that the irreducible modules may occur in V{kXi) (8) ^(Ai) are 
V{{k + l)Ai), V{{k - l)Ai + A2) and V{{k - l)Ai). The multiplicities of them are all less 
than or equal to 1. According to the dimension formula of finite-dimensional irreducible 
modules for a complex finite-dimensional simple Lie algebra (e.g., cf. Page 139 in [H]), 
we have 

diml^((A;- l)Ai + A2) 
n((fc-l)Ai + A2 + 5, a) 

_ gyp 

m rn 

k{k + 2m-2) Yl{k + 2m - 1 - s){s - 1) J] (2m - 1 - j)(j - 1) 

_ s=3 j=3 

m m 

(2m-3) n(2m-l-s)(s-l) n(2m-2-j)(j-2) 

s=3 j=3 

2k{k + 2m - 2){k + m - l){k + 2m - 4)! 



(2m-3)!(A'; + l)! 
k{k + n- 2){2k + n - 2){k + n - 4)\ 
(n-3)!(A;+ 1)! ' 



(2.50) 



where 6 = J2r=i "^r- Moreover, 

dim{V{kXi) O V{Xi)) - dimV{{k + l)Ai) - dimV((k - l)Ai) 
fk + n-3\ fk + n-2\ fk + l+n-3 
= ^ J+^i n-2 )-[ n-2 

'k+l + n-2\_/'k-l + n-3\_fk-l + n- 2 
n-2 ) ~ \ n-2 ) ~ \ n-2 
_ k{k + n-2){2k + n-2){k + n-4)\ 
(n-3)!(A; + 1)! 

= diml/((A;- l)Ai + A2). (2.51) 

Hence all the three modules occur in H'f. ®c Vc- 

Next we consider the case that n = 2m + 1 is odd with m > 3. The irreducible 
modules that may occur in V{kXi) ^^(Ai) are V{{k + l)Ai), V{kXi), V{{k - l)Ai + A2) 
and V{{k — l)Ai). Their multiplicities are all less than or equal to 1 by Lemma 2.2, (2.12) 
and (2.40). Note that V{{k + l)Xi) must occur. Since 

dimy((A;- l)Ai + A2) 
n((A;-l)Ai + A2 + (5, a) 



n «) 

m m 

k{k + 2m - 1)(A; + m - \) Y[{k + 2m - s){k + s - I) Y[{2m - - 1) 

_ s=3 j=3 

mm 

(2m - 2)(m - |) n(2m - s)(s - 1) n(2m - I - - 2) 

s=3 j=3 

_ 2k{k + 2m-l){k + m-l){k + 2m-3)\ 
~ (2m-2)!(A; + l)! 

k{k + n-2){2k + n-2){k + n-4:y. 
(n-3)!(A; + l)! 

= dim(T/(A;Ai) ® y(Ai)) -dimT/((A; + l)Ai) -dimy((A; - l)Ai), (2.52) 

atmosttwoof {\/(A;Ai), V{{k-l)Xi + X2), 1/((A;-l)Ai)} occur. So (2.45) with n = 2m + 1 
follows from the fact 

dimy((A; - l)Ai + A2) > dimy(A;Ai) > dimy((A; - l)Ai). (2.53) 

ii) It follows by a similar argument as that of n = 2m in i). 

iii) It is obtained by a similar argument as that of n = 2m +1 in i). 

iv) The irreducible modules may occur in (8> Vc are the four ones in the right side 
of (2.48) by lemma 2.2. The conclusion follows from the fact 

dimV{{k + l)(Ai + A2)) = (A; + 2)^ dimV{{k - l)(Ai + A2)) = k\ (2.54) 



dimy((A; + l)(Ai) + (A;-l)A2)) 
= dimV((k - l)(Ai) + {k + 1)A2)) = k{k + 2) 



(2.55) 



and 

dim V{k(Xi + A2)) ® ^(Ai + A2) = 4(A; + 1)^. (2.56) 

v) The set of weights occurring in V{2Xi) of o'(3,C) is n(Ai) = {2Ai,0, -2Ai}. The 
statement holds because dim 1^(2 A; Ai) — 2k + 1 for any k > 1. □ 

For convenience, we treat Kr and ^ variables. With these notations, we may write 

n n 
{Er,s - Es,r) ( fj'^j) = {Xrd^cs " ^sd^r + <?r5^, - ?s«9,J /j<?j) (2.57) 

for Er^s — Es^r e o{n, R). 

We define the bar operation on Ac by 



oo n oo n 



y ] y ] '^ri,...,r„,jXl • • • Xn'ij — ^ ] ^ ^ '^ri,---,r„,i^l " " " ^n^jj 0'ri,...,rn,j ^ (2.58) 
ri,-,r„ j=l ri,...,r„ j=l 

For a transformation T on Ac, we define its conjugate operator by 



T{f) = T{f) for / e ^c. (2.59) 

For expository convenience, we shift the sub-indices by —1 when n — 2m + 1 is odd; 
for instance, 

Xt > Xr— 1, Ur ^ Vr—li ^ 'ir—ii f^r ^ '^i — 1) -^r,s ^ -^'r— l,s— 1 (2.60) 

for 1 < r, s < 2m + 1. 

In this way, we always have 

yr = ^ , ym+r- = ^ (2-61) 

and 

for 1 < r < m by (2.25) and (2.30) in the both cases of n = 2m and n — 2m+l. Moreover, 

Ho — ^ ^ Xq and Kq — — if tt- = 2m + 1. (2.63) 

Note that 

y^+r = Wr and Km+r- = for 1 < r < m, (2-64) 

also 

yo — —iyo and = —i^o if n = 2m + 1. (2.65) 



Furthermore, (2.61) yields 

9..- ^"" '7^-+- , 9,^^^ = (2.66) 

and (2.63) says 



9xo = ^9,„. (2.67) 



Thus the Laplacian operator 



(i E if ^ = 2m; 

7^ n^ (2.68) 
^(H'o + E if n = 2m + l. 

In particular, (2.38) holds. Expressions (2.64)-(2.67) also imply 

dys ^9y^: d^^=df^ for 1 < s < n. (2.69) 

A complex module W of o(n, M) is of real type if there exists a real module Wq of 
o(n, ]R) satisfying 1^ = Wq C. In this case, we call Wq the real form of 1^. 

Because we have to deal with the real module finally, we give a lemma to describe 
how to find the real form of a module in H'j^Vc- Recall that we have chosen a Cartan 
subalgebra and positive (negative) root vectors for o'(n, C) in (2.6)-(2.9) and (2.14)-(2.16) 
with sub-indices shifted by —1. A highest (lowest) weight vector of an o'{n, C)-module is 
a weight vector nullified by positive (negative) root vectors. 

Lemma 2.4 Suppose that W C 'H'j^Vc is an irreducible submodule of o'{n,C). Then 
Wq = {h G W \ h = h}, which is an irreducible real o{n,'M.) -module, is the real form of 
W. 

Proof. Recall the action of o'{n, C) on by pi given in (2.32). Observe that 



— —{Vrdys ~ Vm+sdym+r + l^r^K^ ~ K'm+sdKrn+r) ^ (2.70) 



Moreover, 



— ~{yo9yr — Um+rdyQ + KQ^Kr ~ ^m+r^Ko) (2.72) 

if n = 2m + 1. Note for a positive root vector A e o'{n, C), its transpose is a negative 
root vector. The above expressions says that 



A^if) = -A{f). (2.73) 



for / G ^c- 

Now let / be a highest weight vector of W with weight A. Then / is a lowest weight 
vector of some submodule of H'^Vc by (2.70)-(2.72). By (2.36), (2.37) and (2.68), = 
Ti'j.. Moreover, (2.64) and (2.65) imply the weight subspaces 

= (K)-,., = (Vc)-,,. (2.74) 

Thus 

(n'M,^ iK),AVc),,^WM~,-iKVcU. (2.75) 

So / is of weight —A. Since the lowest weight of W is also —A by the highest weights 
listed in Lemma 2.3, the o'{n, C)-siibmodule generated by / must be isomorphic to W. 
But the multiplicity of any o'{n, C) -irreducible submodule in Ti'^Vc is 1 by Lemma 2.3. 
Therefore, f G W. 

Recall that W is spanned by the elements of the form: 

9 = f^lfT,---f^:f (2.76) 

where fr^ is the operator of the r^th negative simple root vector in o'(n, C), G N and 
1 < Tj < m for n = 2m or n = 2m + 1. Then 

m m 

y . — / , (Xt ^ 

(-1)^=1 g = (-1)^=1 A'frJ---f?zf 

rn 

= e«;e«^-.e«:7, (2.77) 

where e,.^, is the operator of the r^th positive simple root vector in o'(n, C). It follows that 
the real and imaginary parts of g are in W whenever g G W . This shows W = Wq + iWo. 
By Lemma 2.1, Wq must be a real irreducible o(n, M)-submodule. □ 

Now we want to find the highest weight vectors of the direct summands of H'/A/c in 
Lemma 2.3. We denote H'^Vc — Wi ® W2 ® W^, where Wj (or Wj) is correspondence 
to the jth direct summand in every equalities in Lemma 2.3 expect the case n — A. Let 
W2 — W2- W2,+ be the sum of the middle two submodules if n = 4. According to 
Lemma 2.4, we have real irreducible o(n, M)-submodules 'Hk,i,'Hk,2,'H'k^^ such that Wj — 
Hkj + i'Hkj for J = 1, 2 and = Ti), 3 + ^7^^ 3. In particular, HkV = Hk,! ® "^^,2 © '^fe.a 
by (2.35) and (2.38). in the case of n = 4, 7Tlfe,2=F similar real forms of VF2,=f- 

Lemma 2.5 i) The following vectors Vj {res. v'^) ofWj {res. WQ and V2,^ ofW2,^ are 
highest weight vectors: 

vi = ViKi = 7^(a;i + ixm+i)''{^i + iwi)> (2-78) 



V2 = -y-2y\~^t^i + y\fi2 for n>A, (2.79) 

V2- = -y2yi"Vi + 2/^«2, V2,+ = -VaVi'^Ki + y^K4 for n = 4, (2.80) 

V2 = -yoyf + «o M ^ = 3, (2.8I) 

m m 

4 = 2(fc - 1) yjym+jyt^Ki - (2/c + n - 4) ^ ^^"^(l/^+jfi:^- + Wj) (2.82) 
for n — 2m, 

m ^ 

4 = -{2k + n- 4)yo?/f-Vo + 2{k - 1)(^ y^y^+j + ^yM'^i^i 

m 

-{2k + n - 4) ^ y'l-\ym+jKj + yjKm+j) (2.83) 

for n = 2m + 1. 

i'ijThe vectors vi, V2 and are solutions of Navier equations (1-6). Moreover, Tik,! 
and 'Hk,2 are subspaces of the solution space. 

Proof, (a) We want to find the highest weight vectors of Wi. The action of positive root 
vectors of o'(2m, C) on are operators: 

yrdy, — ym+sdym+r + l^rdus ~ l^m+sdnm+r (2.84) 

and 

yrdy^^^ — ysdym+r + l^rdnm+s ~ l^sdnm+r (2.85) 

for 1 < r < s < m (c.f. (2.14) and (2.57)). They annihilate v-i. in the case of n = 2m+ 1, 
we have additional operators of positive roots: 

yrdya - yodym+r + K'rdno " Ko^^m+r- (2.86) 

Then (2.84)-(2.86) annihilate Vi. 

Since is harmonic, f i is a singular vector of 7i'^,Vc. The weight of is kXi if n > 4, 
is A;(Ai + A2) if n = 4 and is 2kXi if n = 3. The weight of ni is Ai if n > 4, is Ai + A2 if 
n = 4 and is 2Ai if n = 3. Then the weight of f 1 is {k + l)Ai if n > 4, is {k + l)(Ai + A2)if 
n = 4 and is 2{k + l)Ai if n = 3, which implies that Vi is the highest weight vector of Wi. 
Furthermore, A(t;i) + 6(V^ ■ V){vi) = 6(V^ ■ V){vi). But 

v(^i) = Id^Ay'i) + '^d,r.M) 

k k 
= ^{xi + ixm+i)''~^ + i^^{xi + ixm+i)''-^ = 0, (2.87) 

which implies that vi is a solution of Navier equations. By the invariant property of 
Navier equations under o{n, M), Hk,! is a solution subspace. 



(b) Now we approach the highest weight vectors in W2. It is obviously that V2 G TY'^Vc- 
If m > 1, then 

{Er,s ~ Elm+s,m+r) {^2) 
= iVrdy, - Vm+sdy^+r + ^r^^, - Km+sd^rn+r)i-y2yi~^t^l + yf«2) = (2.88) 

for 2 < r < s < m, 

(-E'1,2 — Em+2,m+l){v2) 
= iyi9y2 - Z/m+29y^+i + - Km+29«^+i)(-y2yf~Vi + y^K2) 

= -y'lKi + y'lKi^O (2.89) 

and 

= (l/r5y^+, - + f^rd^m+s " J (-?/2|/i "^/^l + yiK2) = 0. (2.90) 

Moreover, if n = 2m + 1 is odd, 

(-^r.O — Eo^rn+r)iv2) 

= (yrdyo - yodym+r + K'rdno " /^O^Wr) (-Z/2Z/1 "^/^l + Z/l/^2) = 0. (2.91) 

For n — ?>,V2 — — yoZ/i^^'^i + yi'^o- The operator of positive root vector of o'(3, C) is 

yidy^ - yody2 + f^Ao - i^odK2 (2.92) 
which also annihilates V2- in the case of n = 4, we have that 

iVidy^ - y2dy., + k,i9k4 - K,2d^,J{v2,+) 
= {yidy, + Kid^,){-y4yt^Ki+y^iK4) = -y'lKi + y'lKi^O (2.93) 

and similarly for V2-- 

Note that the weight of V2 is (A; + l)Ai — ai {2(k + l)Ai — ai if n = 3). So it is the 
highest weight of W2. Moreover, the weight of vector V2- (res. ^2,+) is (A;+l)(Ai+A2) — cci 
(res. (A; + l)(Ai + A2) — a2). Hence it is the highest weight of W2- (res. W2,+). 

In the case of n 7^ 3, 

2'=+^V(v2) = 5^i(-(x2 + lXm+2){Xl + ix^+i)''-^) + d^^{Xi + tXm+l)'' 

+dxm+i{-i{x2 + iXm+2){xi + ^) + dxm+2{i{xi + ixm+i)'') 

= -{k- 1){X2 + iXm+2){xi + iXm+lY~'^ 

-i^{k - 1){X2 + ixm+2)ixi + ixm+i)'''^ = (2.94) 

and it is easy to verify V{v2-) = 0. Moreover, 

2'=+V(t;2,+) = dx,(-(ix2 + X4)(xi+ix3)''~^) + dx2(i(xi + ix3)'') 

+dx,{-i{ix2 + X4){xi + ixs)''^^) + dx,{{xi + tx^)^) = 0. (2.95) 



in the case of n = 3, 

+d^,(-(l + i)ixo(xi + ix2)''-^)^0. (2.96) 

Since y'^, y2yi^^ and yoy\^^ are harmonic, A(f2) + 6(V^ ■ V)(f2) = 6(V^ ■ V){v2) = 0. 
Similarly, the same equation holds for f2,^. That is, V2 and V2,zf are solutions of Navier 
equations. Hence 'Hk,2 is a solution subspace. 

(c) We now deal with the highest weight vectors of W^. 

in the case of n = 2m, observe that Em-i,2m — -E'm,2m-i annihilate v'^ and 

= 2(fc - Vjyiym+i+iyi'^'^i - i2k + n- 4:)yiy'['^ Km+i+i - 2{k - l)ym+i+iyiyi'^ k,i 

+ {2k + n- 4:)ym+i+iyi~^Ki - {2k + n - 4:)ym+i+iyi~^ ki + {2k + n - 4)yiy'l-^ k^+i+i 
= (2.97) 

for Z = 1, . . . , m — 1. On the other hand, 

m 

A{2{k yjym+jVt^ -{2k + n- 4)y^+iyiy^2) 

.7 = 1 

= i{2{k -l){k- 1)) -{2k + n-4){k-l) + 2{k - l)(m - l))|/t"^ = 0, (2.98) 



and Ayjy^^'^ = for j m + 1 by (2.68). That is, v'^ e 'H'f.Vc is harmonic , and so it is a 
singular vector of H'^Vc by (a), (b) and Lemma 2.3. 
in the case of n = 2m + 1, we similarly have that 



and 



{yidyi+^ - Vm+i+idy^+i + i^idni+, - Km+Z+l5«^+J(4) = 0. (2.99) 

= -{2k + n- 4:)yryt^Ko + 2{k - l)yryoyi~^i^i - 2{k - l)y^yoyi~^«i 

+{2k + n- A)yoy'[~^Kr -{2k + n- 4)yoyi~^«r + {2k + n- 4)y^yf-^Ko 
= 0. (2.100) 

In addition, 

m ^ 

A{2{k - l){Y,yjym+j + 2?/o)z/i"' -{2k + n- 4)y^+iyiyt^) 

i=i 

= -(A; -l)(2n + 4(A;-2))yi^-2-i(2A; + n-4)(A;- 1)^^-2 = (2.101) 

and Ayjy'l''^ = for j 7^ m + 1 by (2.68). That is, v'^ e "^fc^c is harmonic, and so it is a 
singular vector of H'/^Vc by (a), (b) and Lemma 2.3. Finally, one get v'^ e by checking 
the weight oi v'^. □ 



Set 

m 

^3 = 2 ^ y^ym+jyt^m if n = 2m (2.102) 

and 

m 

4' = (2 X! + yo)yi"'«i if n = 2m + 1. (2.103) 

Since [x^^j;^ — Xsdx^,Yl]^=i^'i] — 0) the irreducible module with v'^ as its highest weight 
vector is isomorphic to V{{k — l)Ai) for n > 4, V{{k — l)(Ai + A2)) for n = 4 and 
V(2(k - l)Ai) for n = 3. Let 

(2fc + n-2)(A; + n-3)(A;- 1) 



2(&-i(2A; + n-4)+A;-l) 



(2.104) 



Lemma 2.6 T/ie vector f 3 = + 0^3 a complex solution of Navier equations. Then 
the complex irreducible o{n, C)-module W3 generated by V3 is a complex solution subspace 
of Navier equations. 

Proof. Suppose n — 2m. We have 



^^3 



= 2{k - 1) X] yjym+jVi - {2k + n - 4)^y'l ^{ym+jHj + Vji^m+j) 
i=i j=i 

n \ ^ 

i=l j=l 

-^—^ — ^^a;,(xi + ix^+i)*^-^^, (2.105) 



^3 



n 1 

: ^ J] a;|(a;i + ixm+i)'''\i - + iXm+i)''~^'^m+i- (2.106) 



Hence 



A(^;) + 6(V^-V)(4) = -^^^-rr^&(^ + n-3)(2^ + n-2)(xi + zx„+i)'=-\i 



2fc 

+ ^^~ ^^ 6(fc + n- 3){2k + n-2){xi+ ix^+i)''-\m+i, 



(2.107) 



A{v;) + b{V^ ■V){v;) = ^{2{2k + n-4) + 2b{k-l)){x, + tx^+, f-\ 

-^(2(2A; + n - 4) + 2b{k - l)){xi + ixm+i)''-\m+i. 



(2.108) 



Assume n = 2m + 1. Then 

2m ^ 2m 

j=0 j=0 

-i + (2.109) 

s=0 



2m ^ 2m 

^3 = ^ + ^Xm+i)''~\ - ^ X] + ^^r^+O^'^^Wi- (2.110) 

j=0 j=0 



Thus 



A(4) + 6(V^-V)(4) = -!l^^6(A; + n-3)(2A; + n-2)(xi + ix^+i)'=-\i 



2fc 

+ ^^~ ^^ 6(A: + n - 3)(2A; + n - 2)(xi + ix^+O'^-'wi, 



(2.111) 



A(4') + 6(V^-V)(4') = ^{2{2k + n-i) + 2b{k-l)){x, + ixm+i)''-\i 

-^(2(2A; + n - 4) + 26(A; - l)){xi + ixm+i)^-\m+i. 

(2.112) 

The conclusion foUows from (2.107), (2.108) and (2.111), (2.112). □ 

Now we can describe the irreducible modules of o(n, M) which are solution spaces of 
Navier equations. 

Lemma 2.7 The linear map determined by 

k+l 

Ipi^ii ■ ■ ■ ^ifc+i) — ^ ^ ' ' ' ^ij-i^ij+i ■ ■ ■ Xik+i'^ij (2.113) 

is an o{n,M.) -module isomorphism from Hk+i to Hk,!- 

Proof. Set f = Xi ■ ■ - x^;^, and then xid^^if) = l2Xi'^^X2~^x''^ ■ ■ ■x'-^. It follows that 

ipixA^if)) 
= ^(/2a:i'^^4'"yi ■■■a^n) 

= kih + l)x['x'i-^x'i ■ ■ ■ <Xi + kik - 1)4'^'4"'4' • • • 4"<?2 

n 

+I2 ljXi~^^X2 • • • x-_lx- Xj^l ■ ■ ■ X^^j. (2.114) 

j=3 



On the other hand, ipif) = IjX^ ■ ■ -x-Zlx- x-^\ ■ ■ ■x^'^^j. Then 

3=1 

n 

{xA^ + <;A^){f) = {xA^ + <iA^){Y^ ljx[' ■ ■ ■ x^zix^'^x^^ ■ ■ ■ x^;:<,j) 

i=i 

= hhx'^x'r^x'i ■ ■ ■ x^-^i + hik - i)x'i+^x'r^x'i ■ ■ ■ x^x2 

n 

-I-/ 7 ^'l + l^'2-l„i3 . . . 'j-l™^J-l . . . J.n^ 

3=3 

= iPixidM). (2.115) 

By symmetry, = A{ilj{f)) for any A G (^/(n, M). Moveover, / is harmonic and 

'^if) = YliQj^ji then all are harmonic. Indeed, we write / = ^ /;(,T2, . . . , a;„)^, 

i=i /=o 
and then / is harmonic if and only if //+2 = —^fi for < Z < A; — 1. It follows 

that Qi = ^ /i+i(x2, . . . , is harmonic. By symmetry, all Qj are harmonic. Hence 

il>{'Hk+i) C Ti-kV is an irreducible module of o(n, M). Thus '4>{Ti-k+i) = 'Hk,i by Lemma 
2.3 and Lemma 2.4. □ 

n n 

Lemma 2.8 i) Denote xj = (k — l) ^ x1<;j — {2k + n — A)xj ^ Xr<ir for 1 < j < n. Then 

r=l r=l 

the linear map ipi determined by 

k-l 

Xii ■ ■ ■ Xii^_^ I >■ ^ ] • • • Xi^_j^Xi^_^_^ ■ ■ ■ Xii^_^Xi^ (2.116) 



is an o{n,M)-module isomorphism from Ti-k-i to 'H]^^- Moreover, — (X] ^r)V' ^■^ o,n 

r=l 

o{n,M^) -module isomorphism from Tik-i to the real form TY^g of the complex o'(n, M)- 
submodule generated by v'^ (cf. (2.102) and (2.103)). 

a) The module is of real type, ip = ipi + Cip2 (cf. (2.104)) is an o{n,W)-module 
isomorphism from Hk-i to the real form (denoted by Hk,3) ofW^ (cf. Lemma 2.6). 

Proof. Firstly, it is clear that 

{xid^^ - x^d^^ + ^i9^2 - ^29^i)(£i) = -^2, (2.117) 

(Xi(9^2 - a^2<9:ri + ?l5f2 - ^2<9a)(x2) = ^1, (2.118) 

{xid^^-X2d^^+<;xd^^-<;2d<i^){xj) = (2.119) 

for j > 2. Moreover, 



99l((Ei,2-E2,l)(C)) = h{ll + l)-Xl + l2il2- 1)^X2 + 



1 ^ ^ "^i 



X2 X2 j—Z •^2Xj 

-hill - i)-^^i - h{i2 + 1)— ^2 - y^Mh~^^3^ 



•^1 "^l X\Xj 



(2.120) 



where C, = x^^ ■ ■ ■ x^^ . On the other hand, 

(£;i,2-£;2,i)(¥'i(C)) 

c c 

= hk—Xi - - 1)^X1 + h — {{Xid,^2 - X2«9a;i + ?l<9a " '?2<9<;i ) l) ) 

a;2 XI Xi 

+12(12 - 1)^^2 - ^1^2—^2 + l2—{{xidx2 - X2dx^ + ^19^2 - q2d^-,){x2)) 



+ y^S^2lj:^Xj - hlj^^Xj) + ^ —{{xidoc2 - X2d^^ + <;id^^ - ?29^i)(5j)) 

(2.121) 



j=Z X2Xj X\Xj .^^ Xj 



by (2.117)-(2.119). Thus (^i((Si,2 - ^2,i)(C)) = (^1,2 - ^2,i)(v^i(C))- By the symmetry 
on sub-indices, ipi is an o(n, R)-module monomorphism from 7ik-i to Ak- If n = 2m, 



<pi{yr) 

k-1 



h[J2[ j{{k-l-r)x'l ^ ''{iXm+iYxi+rix'l ^ ""{iXm+iY ^Xm+i) 

r=0 V ^ / 

hl^[ ^)Xi'^''{iXm+lY{xi + iXm+l) 

r=0 V ^ / 



2 

= -2i(A;-l)v3 (2.122) 

(cf. (2.82)). So </?i is an o(n, M)-module isomorphism from to TtI). 3. It holds similarly 
when n — 2m + 1. 

Lemma 2.8 implies that (p2 is an o(n, R)-module monomorphism from Tik-i to Ak- If 

n 

n = 2m, we denote ^ = ^ a;^, and then 

r=0 V ^ / 
m 

= -4i(A; - 1) J] = -2i(k - l)v'^ (2.123) 

n 

(cf. (2.102)). So (^2 = (5^ 2;^)'0 is an o(n, M)-module isomorphism from Hk-i to the real 

r=l 

form TYfe 3. It holds similarly when n — 2m + 1. Since (fiCH'^^i) — Ws and </?(/) = </?(/) 
for / e we know that is of real type by Lemma 2.4. Moreover, (p{Ti.k-i) is the 

real form of W3. That is (piHk-i) = Ti^.a- □ 



Lemma 2.9 The subspace 

n n 

nk,2 = {J^fr^r I /. e Hk, J^a;,/, = 0}. (2.124) 

r=l r=l 

Proof. Denote by V the right side of above equahty. Note that the hnear map u : 
Ylr=i fr^r ' — ^ Sr=i -^rfr IS an o'{n, C)-module homomorphism from Ac to by Lemma 
2.1. In particular, 

u{Kr) = Vr for 1 < r < n. (2.125) 

Moreover, z/(Re(t;i)) = 'Re{y{vi)) = Re{yi~^^) e Ti-k+i by (2.78) and the classical harmonic 
analysis. Lemma 2.1 also tells us that z/^ is an o{n, ]R)-module homomorphism from A to 
A. Recall that Hk,! is an irreducible o(n, M)-submodule generated by Re(t'i) and 7ik+i is 
an irreducible o(n, ]R)-submodule. So J^l-fy^,^ is an o(n, M)-module isomorphism from Hk,! 
to Hk+i- 

According to (2.61) and (2.63), 

m n 

4^^yrym+r = i^xl, 2y^ = ixl- (2.126) 

r=l s=l 

Thus 

u{v',) = ^~^~\ xl + ■■■ + xlXx, + zx^+.f-' (2.127) 
by (2.82) if n = 2m. When n — 2m + 1, we shift index r i— > r + 1 and have 

^{O = ~lk~ (^1 + • • • + + ix^+2f-^ (2.128) 

by (2.83). Thus v{Re{v'^)) = Re(i/(T;^)) e (x? H h x^)7ifc_i. By the o(n,R)- 

irreducibility of 7Y'j, 3 and {x\ + • • • + x^)Tik-i, ^\h'i_^ is an o(n, R)-module isomorphism 

from 3 to (x^ H h a:;^)'^fc-i. 

By (2.79)-(2.81) and (2.125), u{v2) = and u{v2,±) = 0. Similarly, we have = 0. 

Since HkV = TY^,! 71:^,2 Tijt.s, ker u\n,v = 'TYfe.a, that is, (2.124) holds. □ 

Corollary 2.10 The subspace 

n n 

nk,l+nk,2 = {J2fr^r I fr G Hk, J^^-^ifr) = 0}- (2.129) 
r=l r=l 

Proo/. Denote by V the set in right side in (2.129). Note {d^^ +idx^_^_^){{xi + ixm-i-i)'') = 0. 
By (2.78), Re(fi) e Vfl^fc,!- Moreover, the map a : J2 fr'^r ~^ Yl ^xXfr) is an o(n,R)- 

r=l r=l 

module homomorphism from A to A. The irreducibility of TYj^^i as an o(n, ]R)-submodule 
implies 7^^,! C V. For / e "^^.,2, = A{Y.^rfr) = '^Y^d^Xfr) by straightforward 

r=l r=l 

calculation. Thus iikfi C y. Now aiiikfi) = o-{HkV) ^ {0}. Since ii-k^z is an irreducible 
o(n, M)-submodule, kera|^^^ = {0}. Hence (2.129) holds. □ 



Denote dr , = x.dr^ — Xrdr. and 



( (^3,4 (^4,2 C?2,3 \ 
V 4,2 C^l,3 C^2,l / 



(2.130) 



Then if n = 4, we have that 
Lemma 2.11 The subspaces 

nk,2±^{fenk,2\ vf = ±{k + i)f}. (2.131) 

Proof. We can verify that V commutes with o(4, M). Moreover, 'D{v2,±) = ±(A; + l)w2,± 
by (2.80). Expression (2.131) holds because V2,± are generators of 'Hk,2± as o(4, M)- 
modules. □ 

Now we can get the main theorem of this section. 

Theorem 2.12 Assume that integer n > 3. Let Hk,i — i/j{Hk+i) (cf. (2.113)) and 
'Hk,3 = ^{Ti-k-i) (cf. Lemma 2.8), which are irreducible o{n,'K)-submodules. Take 'Hk,2 
in (2.124), which is an irreducible o{n,M)-submodule if n ^ 4. When n — A, 7ik,2 — 
'Hk,2+®'Hk,2- Ojndi-ik^2± o^re irreducible o(4, R)-submodules characterized by (2.131). Then 
the subspace of homogeneous polynomial solutions with degree k of Navier equation is 
Tik = Hk,! © 'Hk,2 © 'Hk,3, and 

Ak^nk®{xl + --- + xl)Ak-2- (2.132) 

Proof. We have Ti-kj C Tik for j = 1,2,3 by Lemmas 2.5 and 2.6. By Xu's method, we 
can calculate dimTi^ = ndimTi^ = dimTi^ j + dim'Hk,2 + dimTYj^a (The details will be 
given later for technical convenience (see Proposition 3.8)). Since 

Hk = HkV mod (x? + • • • + xl)nk-2V (2.133) 

by the way of our taking subspaces 'Hk,i,'Hk,2 and 'Hk,3 in the theorem, (2.132) follows 
from the facts Ak = Ti-k + {x1 + ■ — h x^)7ife-2, Ak — AkV and induction on k. □ 



3 Bases 

In this section, we will construct some bases of the subspaces Hk,i, 'Hk,2-, 'Hkfi and 'Hk,2^ 
defined in section 2. 

Since 7ik+i = ^ik,! (cf. Lemma 2.7) and Tik-i = 'Hk,3 (cf. Lemma 2.8), we can obtain 
the bases of iik,i and iik,3 by a basis of Hk introduced in [XI]: 

n 

{w(e, /2, . . . , i„) I e e {0, 1}; /„ G N, e + = k} (3.1) 

J=2 



with 

n 

w{e,k,...,Q^ J2 n ^1 (3-2) 

Take e e {0, 1}, ^2, In such that e + Yl%2 Ij — k + 1, and define 

n 

f(e, h,..., In) = il^(w(e, h,..., In)) = fj{e, h, ■ ■ ■ , Q^ij- (3-3) 

Then straightforward calculation shows that 
/i(e, I2, ■ ■ ■ , In) 

n 

= E (^ + 2E^.) H ^1 ^'^^ n^'"" (3-4) 

J=2 



and for j' = 2, . . . , n, 

/j(e, I2, ■ ■ ■ , In) 



j=2 

Proposition 3.1 The set 

n 

{/(e, ^2, ■ ■ ■ , ln)\e ^0 orl; Ij eN; e + J^lj^k + l} (3.6) 

J=2 

forms a basis offCk,!, where the components of f{e, I2, ■ ■ ■ ■, In) o-f^ given by (3.4) and (3.5). 
Similarly, we define 

n 

w{e,k,...,Q^ ^ , X, ' (3-7) 

r.,...,r.^o (i+2eEr,o(tt:;2;:^) 

i=2 

where e e {0, 1}, I2, G N and e + Ej=2 h ^ ^ ~ 1- Thus w e Denote 

n 

g{e,l2,...,ln) ^(p{w{e,l2,...,ln)) = Ef:'(^'^2,---,^n)^-j, (3.8) 

3=i 



then we get that 



gi{e,l2, . . . X) ^ 2^ {{e + 2^rj] 

n 

V ((1. ,^ , (2A; + n-2)(A; + n^3)(A:~l) ,,^ -+^g/--^-,^ 
^ ^^^-^^+ 2(6-i(2A; + n-4) + fe-l) ^%">^ Jl 



Is 27's 



e+2 X; r-j+1 " 

(A;-l)(2A; + n-4)xi H^^''''') (3.9) 

s=2 



and for J = 2, ... , n, 



^j(e, /2, . . . , /n) = J] n — iih - 2?^j) 

J=2 



X ff;t - 1) + + ^ - 2)(fc + n - 3)(fc - ^2'^' h'^-^fr ^i.-.r., 

- (A;- l)(2A; + n-4)xi x^Ylx';'^"^^). (3.10) 

s=2 

Proposition 3.2 The set 

n 

{g(e,l2,...,ln)\e^0 orl, IjEN, e + J^^j^k-l} (3.11) 

j=2 

forms a basis ofHk^s, where the components ofg{e, I2, ■ ■ ■ , In) o-f^ given by (3.9) and (3.10). 

Now we will find a basis of 'Hk,2 by solving the equation in (2.124). Let f = ^ fj^j G 
'Hk,2- Then fj e Hk ior j ^ 1, . . . ,n and /„ = -x~^ J2 ^jfj- We write 

/. = E/n-i,^(^i,---,^n-i)^ for j = l,...,n, (3.12) 



i=0 



where f^_ii{xi, . . . , x„_i) are homogeneous polynomials with degree k — 1. Denote 

s 

^^ = E^l for s = l,...,n-l. (3.13) 

p=i 

Then fj is determined by fn-ifi /n-1,1 "^i^ 

= -A„_i(/^,,,) for < Z < - 2. (3.14) 



Note 



n-l 



n— 1 k 



n-l 



j=\ 1=0 



— —X 



n-l k-1 nj I 



i=l /=0 



Z + 1 U 



(3.15) 



Thus /„ e Tik if and only if 



n— 1 



(a): X]^J-^"-i.o = 

n— 1 rj 
/L\ ^jJn-l,l+3 



n-l 



' l+l 



(3.16) 



) for Z > 



3=1 j=l 

by (3.14) and (3.15). To write down a basis of 'Hk,2, it is sufficient to solve (3.16). 
Lemma 3.3 The following system is equivalent to (3.16). 

n-l 



(a): Xl^^-^^-i-o = 



n—l ra— 1 

W ^ Y^^^AfUl) = -An-l(E^./n-l,l) 



Proof. Note that 

n-l 

-(Z + l)J2^jfU,+s = ii + 3)A„_i(^x,/^,,^J 



(3.17) 



n-l 



n—l 



n—l 



^ (/ + l)5^a;,A„_i(/^_,_,^J = (/ + 3)A„_i(5]x,/^_,,,^J 



n-l 



n-l 



n-l 



n—l n—l 

^ (/ + l)$]a.,(/^i,,+J = -A„_i(5^a;,/^,,,^J. (3.18) 
Set / = in the above equalities, we get that 



n-l 



n-l 



(3.19) 



Thus (3.16) implies (3.17). But 

n— 1 n— 1 

+ 1) = -An-l(E^./n-M+l) 

n—1 n—1 

^ -{l + l)A^.^{J2^.^{fi_,^^_,)) = -A„_i( J] a;,A„_,(/^_,,,_ J) 

i=i i=i 

n—1 ra— 1 rt— 1 

j=l j=l j=l 

n—1 n—1 

^ {I- l)A^.,{J2d.^{fi_,,_,)) = -ALi(E^./n-i,.-i)- (3-20) 

Then by induction, one gets that (3.17(a)) implies (3.16(b)) if / is odd, and (3.17(b)) 
implies (3.16(b)) when I is even. Hence (3.17) is equivalent to (3.16). □ 

Using (3.17) and the condition fj e Hk for j' = 1, . . . , n — 1, we can obtain a basis of 
'Hk,2- For convenience, we classify these base vectors into three disjoint subsets satisfying 
the following conditions, respectively: Let fj be the first nonzero component of the base 

vector / = X) Ml e '^ik,2- 

;=i 

Condition (*): 1 < j < n — 2, and the powers of x„ in fj are even. 
Condition (**): 1 < j < n — 2, and the powers of Xn in fj are odd. 
Condition (***): j = n - 1. 

To find the base vectors satisfying the condition (*), we set /^_]^ ^ = in (3.17(b)) for 
I — 1, . . . ,n — 1. Moreover, for convenience, we write 

k I 
fUo = J2fn-2M^---^^n-2)^ for J = 1 , . . . , n - 1 . (3.21) 

1=0 

Continues the process. In general, we write 

^ I 

1=0 

Then by induction, the following system is equivalent to (3.17(a)) 

/i,o = 



flo-J2fs-iM^---^^s-i)f for s = l,...,n-l. (3.22) 



/^;o = -^7'I?^'-/Io for j>2. 

r=l 



(3.23) 



Now we can write down those base vectors satisfying Condition (*) by (3.23). 

Prom now on, the notations are always nonnegative integers. Assume that fj is the 

first nonzero component of the base vector / = ^ /^Q G ?Yfc,2- Then /jg = by (3.23). 

1=1 

We take fn-ifi be the following monomials: 

/n-1,0 ~ ^1^ ■ ■ ■ ^n-1 ' (3.24) 



where ri + • • ■ + r„_i = k and r^+i + ■ ■ • + r„_i > 0. Moreover, we set 

fLifi-flo forj<g<n-l. (3.25) 

Thus 

n—l n—1 

fn-1,0 = ~'^nH— •+rq,fc(^ ^r9,s)-^ JJ^ (3.26) 



and so 



^=1 ^1 s=l 



2^1 n—l / 21 n—l 



1=0 s=l y(20--. 

for j' < g < n — 1. Then we obtain that: 



(I). The following vectors are the base vectors satisfying Condition (*): 



n—l 



f = J2M'^l-^n'^l^n) (3.28) 
l=j 



for some j G {1, ...,n — 2}, where 



/, = E(-l)'AL.(^p-n^?j (3.2.) 

and /j+i, are given in (3.27) for any nonnegative integers ri, ...,r„_i with ri + 
h r„_i = k and r^+i H h r„_i > 0. 

To get the base vectors satisfying condition (**) and (***), we set fn-i^ = in 
(3.17(a)), and simplify (3.17(b)). Note that it can be written as 

n-2 n-2 

d.^-Af:=ll) + ^n-l{Xn-lf:ill) ^-Y.^M-l,l) - K-lC£^lfn-^,l) (3.30) 

1=1 1=1 

We write 

fc-i 



/n-i,i = E^^'(^i'---'^"-2)ir forl<s<n-l. (3.31) 



x: 

,g!{Xu...,Xn-2] 

1=0 

Substituting (3.31) to (3.30), we get 



left = a.„_,(|:^r%i) + (A„_2+c_j(x._ixi^?r-'4r) 

1=0 ' 1=0 

- i:a+3)^r+Y%+i:^A„_.(^ri^^) (3.32) 

1=0 ' 1=0 



and 



n-2 ; n-2 fe-1 i 

right = -E^-(E^?^)-(^-2+^l-i)(EE^^^^''^ 



n ' 

s=l i=0 s=l i=0 

fc— 1 n— 2 ; fc— 3 2 ^ 

= -EE(^-2^^+^-)(^^^)-EE^^^^2^- (3.33) 

Thus we obtain 

, ^ n— 2 ^ n— 2 

s=l s=l 

for < / < A; — 2, where gfl^^ = = 0. If we define (—1)!! = 0, then 

Z-^Z^y ^) (2s-l)!!(2Z + 2)!!^ + ^S^n-2<^a;J1.5'2s-lJ 

+(-l)'||^AL2(2^o"-^ + E^^^^) - ^ E^^^^Vi (3.35) 



and 



^2"i = E E(-l)"^ |2l)!!(2! + + l)At^a.J(,L) 



n-2 



2/ + 3 



-j^ 7t — ^ 

— ^E^'-^2/+2 (3.36) 



for Z > 0. The above two equahties tell us that is determined by whenever r < n— 2 
and gig"^ under the condition /^_i^o = for Z = 1, . . . , n — 1. Moreover, is determined 
by g'^(xi, . . . , Xn-T) which can be any homogenous polynomial with degree k — \ — s. These 
help us to write down those base vectors satisfying condition (**) and (***). 
Setting g^'^ — 0, similarly as (I), we get that 

(II). The following vectors are the base vectors satisfying Condition (**): 

/ = /i^- + /n-l^n-l - X~'^{Xjfj + Xn-lfn-l)'in (3.37) 

for some j e {1, n — 2}, where 

[|1 / 21+1 \ 

/^=i:(-i)'<-.((lTT)!n<-) (3-38) 



Ill fc-1 r 1 / n-2 



+ Z.^ (2^ + 2)!!(2g)!(2Z + 1)! ^-1^-2 y^^^n-i U ■ 

(-ir^(r„_i-l)!!(r„_i + 2g)!!r„_i! , ^ / r^^+V^'+iTTr 

r„-i!!(rn-i + 2g + 3)!!(r„_i + 2q + 1)1(2/ + i),^^n-i^n-2 ^^^.'^n-i 



n-l 



+ (r„_i + l)r,A^_,A^, (^xt^x^Y^^^ E j )]> (3-39) 

for any nonnegative integers ri, r„_i with ri + • • • + r„_i = A; — 1. 

We set = for r < n — 2 and 

gr'^xY-.-x^-i, (3.40) 
where ri + • • • + r„_2 = k — 1. Then we get 

(III). The following vectors are the base vectors satisfying Condition (***): 

/ = fn-l{<in-l - X'^Xn-l'in-l), (3.41) 

where 

for any nonnegative integers ri, r„_2 with ri + • • • + r„_2 = k — 1. 

Proposition 3.4 T/ie se^ o/^/ie vectors f given in (3.27)-(3.29), (3.37)-(3.39) and (3.41)- 
(3.42) forms a basis of'Hk,2- 

Now we give bases of Hk,2^ in the case of n = 4. Recall that the representation 
operators of negative simple root vectors of o'(4, C) are 

/ai = y2dy^ - ysdy^ + K2d^^ - Ksd^^ (3.43) 

and 

/aa = Vidy, - y^dy^ + n^d^^ - n^d^^. (3.44) 

Then 

7Y.,2- = span {Re(/^j;^(t;2,-)), Im(/^ (^2,-)) | r, . e M} (3.45) 
by Lemma 2.4. Observe that 

/aiK-) = {y2dy, + «29«J"(-y2Z/l"Vi + yt«2) 

= (z/29; + rz/2^-'«29;-^9,J(-y2yf-'«i + Z/f«2) 

= (H ^ - 1 - j)(-Z/2+'l/i-'-'''^i + Z/2Z/t-''/^2) (3.46) 

3=0 



for r < k — 1, and 

fa,i-y2^Vi-'-''^i + f2yt''^2) 
I -{y^dy, - y^dy,)Y2''y\-'"' - ^s{yA^y, - y^dy^f-'ylyi-'' ' 

{yidy, - y3dy^yy^2yi \- is{yAdy^ - y^dy^y-'^yl'^^y^ ^ 



1 

2 



-l.,r-,fc— r 



i^y^dy.-y^dy^yyl-^^yl ^ - s(y45j;, - ys^j^J" ^l/ai/i 
V ^(|/4a,, - |/3a,J%2^|/r^' - s{y,dy, - y^dy^y-'yl-^'yl-'-' 

for < s < A; + 1. Denote 



^+(r, 5) = 2'^Re((|/4a,, - y^dy.yylyl-'-) 



and 



Then the real part of (3.47) is 

/ -g+{r + 1, s) + sg_{r, s-l)\ 
1 5f+(r,s) + s5f_(r + l,s-l) 
2fc+i ^_(r+l,s)-sc/+(r,s-l) 

\ -g-{r, s) - sg+(r + 1, s - 1) J 



v{r, s) — 



and the imaginary part of (3.47) is 



2k+l 



I -g- (r + 1, s) - sg+{r, s - I) \ 
g-{r,s) - sg+{r + l,s - 1) 
-t/+(r + l,s) - sg^{r,s- 1) 
\ g+{r,s)- sg^{r + l,s-l) J 



Moreover, a straightforward calculation shows that 
i)IfO<r + s<A; and < r < s, then 

1=0 p=0 <j=0 

-/ fc— ) — s—p 



ii) If < r + s < A: and < 5 < r, then 

s k—r—s r—s 



■+p+q 



1=0 p=0 q=0 
„2 , J2\l(J2 



^k — 
.s - I . 



x(^x'^ -\- X^) (^2 + X4)* ^ ^ * 



iii) U r + s > k and < r < s, then 

r+s—k s—r 



s r+s—K s—r \ 
Z=0 p=0 g=0 ^ ^ 



x(a;2 + x^)^^-^-^+'^-2 



^2 I ™2Nr-Z r+s-fe-p g p s-r-g 
Ji"2^ ~r 4/ 1 Ju'2^Ju^Ju^ 



I 



^k — r — s 
K P 



^k — r — s 



''r + s — k 
K P 



(3.47) 



(3.48) 
(3.49) 



(3.50) 



(3.51) 



s — r 



(3.52) 



r — s 



(3.53) 



s — r 



(3.54) 



iv) If r + s > A; and < s < r, then 

s r+s—k r—s 



1=0 p=0 q=0 ^ ^ ^ 



r + s — k\ — s 



Now we have the following base vectors of 'Hk,2-' 
i) 



r ^ .f I < r < If] / III + 1 < r < A; - 1 



f ^ -f J < r < [f 1 / III + 1 < r < A; - 1 

/ = ^(r, s), If ^ < ^ < jliij or I ^ ^ , (3.57) 



ii) 



iii) 

/ = I '\ ^! 'I ^ ° for r = . = ^ and is odd. (3.58) 

Then by (3.45), we get that 

Proposition 3.5 The set of vectors f given in (3.56)-(3.58) forms a basis o/7ife,2-- 

Remcirk 3.6 We observe that there exists a linear isomorphism between vector space 
Hk,2- andHk,2+ 

a : 7l!fc,2- — ^fe,2+, where a{x]) = x\, a{x-2) = X4, (7{x-^ = X3 anc/ (7(a;4) = X2. 

These give a basis ofi-ik,2+ by easily interchanging X2 and X4 in a basis ofi-ik,2-- 

At the end of this section, we will use Xu's method in [XI] to construct a uniform 
basis of the polynomial solution space of Navier equations, whose the cardinality was 
pre-used in the proof of Theorem 2.12. It is different from those bases given above. It is 
not listed in accordance to the irreducible summands of the polynomial solution space. 
Xu's method is also critical to solve the initial value problems of Navier equations and 
Lame equations in the next section. 

Lemma 3.7 [Xl] Suppose that is a free module of a subalgebra SS generated by a 

00 

filtrated subspace V = [J Vr (i.e. Vr G Vr+i). Let To be a linear operator on with 

r=0 

right inverse Tq" such that Tq{^) C Tq{^) C ^ and To{rji'r]2) = T'o(?7i)?72 for rji G 
r]2 e V. Let Ti, . . . ,Tm be linear operators on such that Tjiy) C V , Tj{fQ — fTj{Q 
/or j = 1, . . . , m, / e C e •2^- If T^{h) = with he ^ and g eV, then 

00 m 

^-Y.^Y.^^orTsy{hg) (3.59) 

i=0 s=l 



is a solution of the equation 



(3.60) 



Suppose TjiYr) C Vr-i for j = 1, • • • , m, r e N, where V^i — 0. Then any polynomial 
solution of (3.60) is a linear combinations of the solutions of the form (3.59). In partic- 
ular, if TrTs = TsTr, ToTj = TjTo and TqTj = TjTq for any j, r,s e {1,..., m}, then u 
can be written as follows: 



u 



= E 



2i H hi, 

^1) ■ ■ ■ ) 



Note that Navier equations (1.6) can be written as the following form 

(r2-rori-r2)(«) = o, 



where 



To = - 



/ n 
6+1 ^ Xj 









Ti = - 



6(9.. 



6+1 X" 






\ 



It 



j=2 



X2 '-'Xz 

n 



bd^codxr 



bd^^d^^ 



\ 



i=2 



bdxA 



j=2 



bdx„dx2 
where 6 = (ti + i2)/ ii- Set ^ = M[a;i]/„ and 

n 
r=l 

By Lemma 3.7, any polynomial solution of (3.62) is a linear combination of 

oo 



m=0 



(3.61) 



(3.62) 



(3.63) 



(3.64) 



(3.65) 



(3.66) 



where h — x\In and g & V . In order to write down these solutions explicitly, we need to 
calculate (Tq~)'"(Ti + TqT2)'^. Recall that £'r,s is the nx n matrix whose (r, s)-th entry 
is 1 and the others are 0. We define the hnear operator J^^^^ on ]R[xi, . . . , x„] by 

1 



J(xi) ri + i 



~n+l„r-2 



(3.67) 



We take 
Then 



To" = / In. 



(3.68) 



-(Ti+T„-T2) 

{{b + 1)-' E [ )E,,, + Y.{^E,,r + hd^^Er,,) 

n II „ n p 



(3.69) 



r=2 j=2 



r,s=2 



Observe that J^^^y 8^2, ■ ■ ■ , d^^ commute pairwise, and so the entries of (Tq )'"(Ti+Tq 



are polynomials in M[j^^^^, dx2, • • • , dx„]. In order to use linear algebra, we replace these 



operators by real numbers as Xu did in [X2]. Indeed, if we let 



B{ai, 0,2, 



' b+1 b+1 b+1 

ba2 baial + ^ 6010203 
60,3 6010302 6aia| + ^ 



y ban baiana2 6010^03 



ban \ 

6aia2a„ 
601030^ 

6aia^ + ^ y 



with ds e M, 



then 



77 = and ^ = oiTy, 

i?(^y" ,a2,...,a,„) = -(Ti + To-r2). 



Observe that for m > 1, 



with 



6+1 



In-l 



ba2 



In-l 



/b+1 



bas 
/b+1 



60102 + ^ 6010203 

= 6040302 6aio| + ^ • 



yb+i^ 
601020^ 

ftaiOsOj, 



V ^^fei ^^i«n^2 baidnds ■■■ bdiaf^ + ^ J 
is symmetrical (note that 6 + 1 = (2ti + t2)/ti > 0). Moreover, we denote 

a = (6 +1)2 + 1 and = (6 + 2)^^^ ^ 4(6 + 1)77. 

It can be proved that the eigenvalues of A are 

^1 = 4?^ and ^2="^"^^ 



2(6+1; 



2(6+1) ' 



(3.70) 



(3.71) 



(3.72) 



(3.73) 



(3.74) 



(3.75) 



(3.76) 



where the multipUcity of ^ is n — 2. Recall that = (0, 0, 1, 0, 0)"^. We can take 
orthonormal eigenvectors Vj = {0,Pj^2, ■ ■ ■ iPj,n)^ j = 1) ~ 2 corresponding to the 
eigenvalue ^, a unit eigenvector 



^/w + {b + 2)CV^ ^ V2r){w +{b + 2)CV^) 
corresponding to the eigenvalue 9i, and a unit eigenvector 

V. = + ± _^M±MZ^,.. (3.78) 

v/w-(& + 2)eV^ ^ V2?7(^-(& + 2)eV^) 

corresponding to the eigenvalue 62. Setting P = {vi, - ■ ■ , Vn) and 

/ a^-2 \ 

J = ^1 , (3.79) 
V ^2 / 

we have that 

A"' = PJ'^P' = P ■ diag(r , • • • , r , OT, 0T) ■ = {Cr,s)nxn, (3.80) 

where 

Cl,j = Cj^l = Pl.n-lPj.n-l^T + Pl,nPj,nOT (3-81) 

for 1 < j < n, and 

n-2 



Dm 
''2 



^ V PrjPsjC Pr,n—lPs,n—10l Pr,nPs,n02 
= +Pr-,n-lP.,n-l(^r " r) + P.,nPs,n(^2"^ - D, (3-82) 

for r, s e {2, . . . , n} by the the fact PP^ — /„. Substituting (3.76) into the above two 
equations, we have that 



■2(b + 1)// + bJ^)"' 2(h + 1)// (a^ - bJ^Y 



ro + (6 + 2)^V^ 2™(6 + l)"^ ro-(6 + 2)^V^ 2'"(6 + l)"» 

^ - + 2) (^^^ a™-^6^C"'+'w^) 

(3.83) 



2"*(6 + l) 



^2(6+1)77 a^((& + 2)^ + v^) (aC + 6V^) 



v/^T(6T2)e^ v/2r7(tz7 + (6 + 2)^V^) 2-(6 + l)- 
V2(& + 1)^ a,((6 + 2)e- V^) (ae-6x/^)' 



^w-{h^2)i^ ^2i]{w-{h + 2)i^) 2-(6 + l)- 
r y ('^\ oT'-'WC-^w'^ (3.84) 



for 2 < r < n, 



-'r,s 



-( 



a^((6 + 2)^ + v^) a,((6 + 2)^ + 



V2r/(w + (6 + 2)^x7^) V2^(w + (6 + 2)Cx/^) 



a,((6 + 2)e- V^) a,((6 + 2)e- V^) 



or 



V2r/(ti7 - (6 + 2)eV^) V'2;K^^T^T2je^ 
a,((6 + 2)^ + v^) a,((6 + 2)e + v^) (a^ + 6^^) 



^277(07 + (6 + 2)ex/^) V2^(tu+(6 + 2)eV^) 2-(6+l)- 
a,((6 + 2)e- v^) a,((6 + 2)e- v^) (a^ - 6v^y 



V27y(tz7-(6 + 2)ey^) v/2r/(ti7-(6 + 2)eV^) 2-(6 + 1)' 



2b- 



2tJ 



for 2 < r, s < n with r ^ s, and 



61(77 -a^)r"' + 



2"*(6 + 1)^ 



+(6 + 2)5^("')a-^Ve 



2b 



for 2 < r < 7^. 

For convenience, we denote 



(3.85) 



(3.86) 



/(m,s) = 



A'{b + iy{b + 2)2'-2^a™-2'-62'- 
2"»(6+ I)'" 



r \ I m 
s) \2r 



(3.87) 



and 



^(m,s) = 2^ 



Thus we have 



(arS-)(ai,a2,...,a„) = 



2"^(6 + 1)' 



sy/ V2r + 1 



Am „ 
"1 '-2,n 



for 77T, > 1, where 



s=0 



(3.88) 



(3.89) 



(3.90) 



^ 2^(m,s)a^-2-^77— ^a,-, 



s=0 



(3.91) 



[^1 

VbTTa'['cj,i^ ^ 2(6 + l)^(m,s)a?'"-2"-V~'"^% (3.92) 

s=0 

for 2 < J < n and 

Ifl 

aTcj,i = 5j,ia'rrj^ + j](-(5,,o + /(m, s) + {b + 2)g{m, s))a'r-'r"%ai (3.93) 

s=0 

for 2<j, l<nhy (3.73), (3.83)-(3.89). According to (3.72), 



Uj{€, I2, ■ ■ ■ , In) 

of Navier equations, where 



\ U]{e,l2, ...,ln) 



n CO U 2 Jl 

ul = + J]^(-l)-(/(m,s)-(6 + 2)5(m,s)) 

q=2 m=l s=0 

„e+2m-2s 



- (.X-2.)! ^"'"(n4-). (3.97) 



9=2 



°° " 2 J e+2m-2s-l 



T m — 1 T 



(3.94) 



Set 

A2,n = EC- (3-95) 

r=2 

By (3.66) and (3.89)- (3.94), we have the foUowing solutions 



(3.96) 



= E + ^)^^ re + L -2s- IV ^^" n (3-98) 

m=l s=0 ^ ^' q=2 

for 2 < r < n, 



00 11 2 i e+2m-2s-l ™ 

= E E (-l)-2,(m, .)/, _ A^r-^(^7^ n 4-) (3-99) 

m=l s=0 ^ q=2 

for 2 < j <n, 

m=l s=0 

e+2m-2s " 

^e+'2^-2.)! ^""-''""-^"n4') (3.100) 

^ ' g=2 



for 2 < r, j < n with r 7^ j, and 



e+2m 



-'in4 + E(-ir77Tw^2'-'n4') 

q=2 m=i le + zm;. 

00 [f] 

+ E E(-l)"^(-'^.,o + f{m, s) + {b + 2)g{m, s))lj{lj - 1) 

m=l s=0 

e+2m-2s " 

g=2 



X 



for 2 < j < n. 



Theorem 3.8 T/ie set { Uj{e, I2, . . . Jn) \ e — or 1; Ir j = 1, . . . , n} forms a basis 
of the space of polynomial solutions for Navier equations. In particular, 

dimHk — n ■ dimHk- (3.102) 

4 Initial Value Problems 

In this section we will deal with the initial value problems of Navier equations and Lame 
equations using Lemma 3.7 and Fourier expansions. 

Recall that the initial value problem of Navier equations is as follows: 

tiAfw) + (ti + t2)(V^ • V)(«) = 

U{0, X2,..., Xn) = go{x2, ■ ■ ■ , Xn) (4-1) 
U:ci{0,X2, ...,Xn)= gi{x2-, ■■■■,Xn) 

with xi e R and ] for r = 2, . . . , n, where 02, . . . , On are positive numbers, 

and the jth component gl of g^ is a continuous function for e = or 1 and j — 1, . . . ,n. 
For convenience, we denote 

kl = ^, k^ = {kl kif for k = {k2, knf e (4.2) 

dip 

and 

n 

X —— (^2? • • • 5 "^n) 5 k^ ' ^ ~~ ^ ^ k^Xrf. (4.3) 

r=2 

By lemma 3.7, we have that 

00 

Uxi. . . . , ^n) = E(^o~)"^m + T^T2r{x\g:} (4.4) 

m=0 

and 

oo 

V^i(xi, . . . ,x„) = E(ro")"(7^i + ro-r2)'"(xia,,(0o(O,a;2, . . .,xn))) (4.5) 

m=0 



are solutions of Navier equations. Denote 

u{xi, . . . ,Xn) = $o{xi, ...,Xn) + 01 (^i, . . . , a;„) - ifi{xi, . . .,Xn). (4.6) 

Then by (3.90)-(3.94) and superposition principle, the function u{xi, . . . is the solu- 
tion of (4.1). Now we give the explicit expression of u{xi, . . . ,Xn)- For convenience, we 
write that 

n n 

l = Y,€,,, ^P, = Y,^{,,. (4.7) 

3=1 j=l 

Take the Fourier expansions of 



Qe^Yl Yl {H0)cos27r{P ■x)+4{k)sm27r{P ■x))qj, (4.8) 



1 



2"-2a2 ■■■a 
1 



where 

biik) 

Substituting (4.8) into (4.4), and using (3.90)- (3.94) again, we get that 

°o e+2m-2s 

- E E(E(-l)^4— (/(m, s)-ib + 2)g(m, .)) ^^^^ _ 



/•aa /-an 

-/ ■■■/ gl{x2, . . . ,Xn)cos2n{k^ • x)dxn- ■ -(1x2, (4.9) 

'^i J —a2 J —an 

/a2 ran _^ 

■■■ gi{x2, ■ ■ ■ ,Xn)sm2TT(k^ ■ x)dxn- ■ ■dx2- (4.10) 
■02 J -an 



n 

x{bl(k) cos27r(^t • x) + cl(k) sin27r(^t • x)){J2iH^)^r''' 

r=2 



" ' 2 ii e+2m-2s-l 

Xi 



Z=2 s=0 ^ ' 

n 

x{b[{k) sm2n{P ■ x) - c[{k) cos 27: {k^ ■ x)){J2{kUfr~'~^), (4-11) 

r=2 

°o I^^ll c+2m-2s-l 

^. = E E(E(-ir^"-(^+ito('».^K (,^2'm-2.-i)! 

n 

x(bl(k) sin27r(^t • x) - cl(k) cos27r(^t • f))( J](fc^7r)')"^-"-^ 

e+2m n 

+4^ ^^ _^^^^, (6^(fc) cos27r(fct • f) + c^,(fc) sin27r(fct • x))(^(A;^7r)=')'" 

+ E E(-l)^4--7r2^J^|(-5,,o + /(m, s) + (6 + 2)^(m, ^)) ^ - 2^)! 

n 

x(6^(A;)cos27r(A?t -f) + 4(A;)sin27r(A?t •f))(5];(A;^7r)2)'"-^-^) (4.12) 



r=2 



for J = 2, ... , n, 

oo n ItI , -..suii-s 07,_7,t ,,.l-2/)i-2,s 

- E E(EE - (^^^M'".-)) (iX-2.)! 

0-<fcsN"~i m=0 1=2 s=0 

n 

x{b\^{k) sm27r(^t • x) - cos27r(P • x))C^{kl7rf)"'-' 

r=2 

n 2m-2s 

+ E E (-1)-+^4—-2KW)^^K^) . 2^^2.11 

1=2 s=0 ^ '' 



x{bl{k) cos27r(^t . x) + cj(fc) sin27r(^t • f))( ^^(^Tr)')"^-^-^), (4.13) 

r-=2 

and 

oo n 2m-2s 



E E(EE(-ir^4— •2^A^^]^K^)7^^^ 

0-<fc£N"~i m=0 1=2 s=0 ^ ' 

n 

x{b\^{k) cos27r(^^ • x) + 4(^) sin27r(^t • f))(^(A;j7r)')'"-^-^ 

r-=2 

™l+2m " 

+4"^ • 267rA;J— rribUk) sin27r(^t • f) - cl(k) cos27r(P • f))(^(A;]:7r)2)' 



" 1 2 Jl l+2m-2s 

+ E E(-l)^4— • 27r3^J(fct)2(_^^^^ + ^(^^ ^) + + 2)5(m, ^)) ^ - 2.)! 

n 

x{bl{k) sin27r(A^t • f) - cl{k) cos27r(A^t • :?))(X]('^^^)')'""'"') (4-14) 

r=2 

for j = 2, . . . ,n. Thus, as we mentioned, by superposition principle and Fourier expan- 
sions, we get 

Theorem 4.1 The solution of (4-1) is 

n 1 

U{xi, ...,Xn)= EE ^'(^1' ■ ■ ■ ' ^n) - 'i/^iixu 2;„))?;, (4.15) 



where (f)[{xi, . . . , Xn) and ip[{xi, . . . , x„) are defined by (4-ll)-(4-M)- The convergence of 
the series (4-15) is guaranteed by Kovalevskaya Theorem on the existence and uniqueness 
of the solution of linear partial differential equations when the given functions in (4-1) are 
analytic. 

The initial value problem of Lame equations is as follows: 

ir«-r^A(ir)-(v^-v)(ir) = o 

{t(0, xi, . . . , Xn) = /io(a^i, ■ ■ ■ , a^n) (4.16) 

Mt(0, Xi,...,Xn) =hl{xu..., Xn) 



with Xr G [—ttr, flr] for r = 1, . . . ,n, where ai, . . . , a„ are positive numbers and the jth 
component of is continuous function for e = or 1 and j = 1, . . . ,n. 
Set 

Ti = dfln, T2 = b-^AIn + H, where i/ = V^-V. (4.17) 
Then by Lemma 3.7, the set 

oo 

{T.('^innT,"^)m l i; e = O or l} (4.18) 

m=0 

spans the polynomial solution space of Lame equations. Note that 

= &-'"(A"^4 + {{b + - l)A'^-^H) (4.19) 

for m > 1 and T2 — In- 

For convenience, we denote 

^ ^ (^t^ ...^ fo, ^ ^ (^^^ ^ e N^ (4.20) 

and 

n 

X = (xi, . . . , a;„)^, • X = a~^A;j.a;r. (4-21) 

r=l 

Take the Fourier expansions of h^: 

n 

he^^ {bi{k) cos27r{P ■x) + 40) sm27r{P ■x))c;j, (4.22) 



where 



VJk) = / ••• / hUxi, . . . ,Xn) cos27r(k^ ■ x)dxn- ■ ■ dxi, (4.23) 

2^* ^ai---anJ-ai J -an 

= 7^— i / ••• / hl{xi,...,Xn)sm27r{P ■x)dxn---dxi. (4.24) 

2" 'ai---anJ^ai J -an 

Note that 

CXD 

0, = ^(rr)-(f(rr)(^,)) (4.25) 

m=0 

are solutions of Lame equations for e = 0, 1. Then by superposition principle, the vector 

1 

u{t, Xi,...,Xn) = Y 0e(i, xi,..., x„) (4.26) 

e=0 

is also a solution. Moreover, one finds easily that 

— * 

u{0, Xi,..., Xn) = hQ{xi, ...,Xn) (4.27) 

and 

— * 

ut{0, xi, . . . , x„) = hi{xi, ...,Xn). (4.28) 
Thus it is the solution of (4.16). Substituting (4.22) into (4.26), we get: 



Theorem 4.2 The solution of (4-16) is 

1 n oo (-,\m4.e+2m 

^ = E E E E 2m)! ^^^'^^^ • + ^-^^^ • 

n n 

x(5^(27^^t)2)n^ + 1)- - l)A^k]k\ 

r=l 1=1 

n 

x(b[(k) cos27r(^t • x) + c[{k) sm27r(^t • x))C^{2T:klf)'^-^)qj (4.29) 

r=l 

T/ie convergence of the series (4 -29) is guaranteed by Kovalevskaya Theorem on the exis- 
tence and uniqueness of the solution of linear partial differential equations when the given 
functions in (4-16) are analytic. 

Set 

r 

5; = (0,...,0,4^---a;l",0,...,0)^ (4.30) 

and 

T^{gr)^{g,,...rgnf. (4.31) 

Then by (4.19), we have 

n n 

~g, = 5,>-"(A™(n4l + iMib + 1)™ - l)A-\x;'xj'l[x[')). (4.32) 

s=l s=l 

Thus by (4.18). we obtain 



Proposition 4.3 The set 

If 

j=l m=l 
n 

xA'^-\x;^xfY[xi^))qj\r^l,...,n, e = or 1, e N} (4.33) 

s=l 



s=l 7 = 1 m,=l ^ ' s=l 



forms a basis of polynomial solution space of Lame equations. 

Acknowledgement: I would hke to thank Professor Xiaoping Xu for his advice and 
suggesting this research topic. 
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